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Influence of disorder on the temperature of superconducting transition (T c ) is studied within the 
cr-model renormalization group framework. Electron-electron interaction in particle-hole and Cooper 
channels is taken into account and assumed to be short-range. Two-dimensional systems in the weak 
localization and antilocalization regime, as well as systems near mobility edge are considered. It is 
shown that in all these regimes the Anderson localization leads to strong enhancement of T c related 
to the multifractal character of wave functions. 
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Soon after the development of the microscopic the- 
ory of superconductivity (SC) by Bardeen, Cooper, and 
Schrieffer (BCS) 0], the question of influence of disorder 
on SC attracted a great deal of attention. It was found 
0, HI that the diffusive motion of electrons does not af- 
fect essentially the temperature T c of superconducting 
transition, i.e the mean free path does not enter the ex- 
pression for T c . This statement is conventionally called 
"Anderson theorem" . 

Effects of disorder-induced Anderson localization Q 
on SC were considered in Refs. 0, Q. It was found that, 
within the BCS approach, the SC in a disordered sys- 
tem persists up to the localization threshold and even in 
the localized regime near the Anderson transition. Fur- 
thermore, Refs. 0, 0] came to the conclusion that the 
mean-field T c in these regimes remains unaffected by dis- 
order (i.e. the Anderson theorem holds). In a parallel line 
of research, it was discovered 0-0 that an interplay of 
long-range (1/r) Coulomb interaction and disorder leads 
to suppression of T c . These ideas were put on the solid 
basis by Finkelstein [1(1 EH who developed the cr-model 
rcnormalization-group (RG) formalism. 

Recently, Feigelman et al. 11211 . motivated by a large 
body of experimental data [13|, returned to the prob- 
lem of interplay of disorder and SC. They found that the 
eigenfunction multifractality near the localization thresh- 
old strongly affects properties of a superconductor. Their 
remarkable finding is that T c is dramatically enhanced: 
its dependence on the coupling constant is no more expo- 
nential (as in the conventional BCS solution) but rather 
of a power-law type. This result was obtained on the 
basis of the BCS-type self-consistency equation, with 
Cooper attraction being the only interaction included. 

In this paper we reconsider the problem of interplay 
of SC and localization in the framework of the cr-model 
RG. We take into account the interaction in all channels 
(singlet and triplet particle- hole, and Cooper). The ke 
assumption that distinguishes our work from Ref. 



is the short-range character of interaction that physi- 
cally corresponds to a strong screening of the long-range 
Coulomb interaction. We consider first two-dimensional 
(2D) systems in the weak-localization (WL) and weak- 
antilocalization (WAL) regimes and find a strong en- 
hancement of T c by disorder. We show that this effect 
can be traced back to multifractality of wave functions. 
We then extend the analysis to the vicinity of Anderson 
transitions in 2D and 3D (three-dimensional) systems. 

We begin by considering a 2D system of the orthog- 
onal symmetry class (i.e. with preserved spin-rotation 
symmetry). In the cr-model formalism the problem is 
characterized by four running couplings: the dimension- 
less resistance t = 2/irg (where g is the conductivity 
of the system in units of e 2 /h) and three interaction 
constants Tj corresponding to singlet particle- hole (IV), 
triplet particle-hole (r t ), and Cooper (r c ) channels [lCj. 
The RG treatment yields a set of five coupled equations 
for these couplings and the field renormalization con- 
stant z. It is convenient to switch from Ti to normal- 
ized coupling constants ji = Ti/z (typically, 7 S ,7 C < 
and 7 t > 0); then the equation for z does not affect the 
remaining equations. The case of long-range Coulomb 
interaction corresponds to j s = — 1. 

Assuming weak short-range interaction, |7,| <C 1, we 
have obtained the following RG system of equations (see 
Supporting Material [3): 
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Here y = InL, where L is the running RG length scale. 
We measure lengths in units of the microscopic scale 
where the RG ([I])-© starts. All the equations are writ- 
ten to the leading order in t, 7, <C 1. The first term on 
the r.h.s. of Eq. ([TJ represents the effects of disorder on 



2 



interaction, while the second term is the Cooper renor- 
malization of j c . The r.h.s. of Eq. <j2j) describes the WL 
effect. We discard the Altshuler-Aronov-type contribu- 
tion in Eq. ([2]) (describing renormalization of disorder by 
interaction) since it is of higher order in 7.; . 

Let us analyze the RG flow governed by Eqs. ([1]) and 
([2]). Equation ([2]) decouples from the rest, yielding 



i" 1 (y) = *o 1 -v, 



(4) 



where the subscript refers to the bare value of the cor- 
responding coupling. This is the usual WL behavior. In 
the absence of Eq. ([I]) it would imply that strong Ander- 
son insulator emerges at the scale y = t$ . We now turn 
to Eq. ([1]) . Let us assume that the disorder is sufficiently 
strong compared to the interaction, so that to 3> |t»,o|- 
Then at the initial portion of the RG flow we can ne- 
glect the Cooper renormalization term, which leaves us 
with a linear system of equations. The corresponding 
3x3 matrix has two eigenvalues: a positive one, A = 2t 
and a doubly degenerate negative one, A' = — t. There- 
fore, as RG starts to operate, the vector formed by three 
couplings quickly (at y ~ 1) approaches the eigenvector 
corresponding to A, i.e. — 7 S = 74 = 7 C = 7. Projecting 
the system (p} onto this eigenvector, we get 



dj/dy = 2t 7 - 27 2 /3. 



(5) 



We have checked that the neglected contributions do not 
affect the results in any essential way [14]. 

We will assume the initial value 70 = (— 7s,o + 37 t ,o + 
27c,o)/6 to be negative which will imply SC (or at least 
a tendency towards it) [l5j |. This is, in particular, the 
case when the dominant bare interaction is the Cooper 
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attraction 7 Cj o < 0. Solving then Eq. ((5j), we find 
that there are two distinct situations. If |7 | <C t^, * ne 
resistance t reaches a value of order unity when the in- 
teraction is still weak. This means that when the scale 
is further increased (i.e. the temperature is lowered), 
the system becomes an insulator. On the other hand, if 
l7o| 3> to, the RG flow develops a superconducting insta- 
bility (blow up of the interaction) at a scale where the 
resistance is still small, t <C 1. This RG scale determines 
the temperature T c of superconducting transition, 



T c ~ exp{-2/t } , 
which is much higher than the clean BCS value 



T. 



BCS 



exp{-l/|7 c ,o|} 



(6) 



(7) 



in the considered regime of sufficiently strong disorder, 
to \ji,o\- When the dimensionless resistance becomes 
smaller than the interaction, Eq. (|6]) crosses over into 
Eq. (|7|). We thus find that T c shows a non-monotonous 
dependence on the disorder strength and gets strongly 
enhanced (by a parametrically large factor in the ex- 
ponential) in the intermediate range of resistivities to, 
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FIG. 1: (Color online) Temperature dependence of re- 
sistivity t near SIT in the 2D orthogonal symmetry 
class from numerical solution of Eqs. JT]) and (J2J) for 
7s0 = -0.005, 7 to = 0.005, 7 c0 = -0.04, and t = 
0.065,0.075,0.085,0.095,0.10,0.105,0.11,0.12 (from bottom 
to top). The point where 7 C diverges determines T c provided 
t < 1. Inset: dependence of T c on the bare resisitivity to- 



see Fig. Q] For given interaction strength, T c is the 
largest when the system approaches the superconductor- 
insulator transition (SIT) that takes place at t\ ~ |7 |. 

It is important to emphasize a relation between this 
RG analysis of the interacting problem and the multi- 
fractality of wave functions in the non-interacting the- 
ory. The wave function (or, equivalently, local density 
of states) multifractality is a hallmark of criticality in- 
duced by Anderson localization [l6| . It implies anoma- 
lous power-law scaling of moments (and, more generally, 
of correlation functions) of wave function amplitudes. In 
the field-theory language the corresponding exponents 
are scaling dimensions of composite operators. In the 
fourth order in wave function amplitudes -0, which corre- 
sponds to the second order with respect to the er-modcl 
field Q ~ ipip there are two such operators. The domi- 
nant one has a negative scaling dimension (i.e. it is RG 
relevant), A2 < 0, which is the most famous representa- 
tive of the family of anomalous dimensions A„ describing 
the wave function multifractality spectrum [16l4l8| . The 
second operator has a positive dimension (12 > 0, which 
means that it is RG-irrelevant. It is worth mentioning 
that, despite the RG-irrelevant character, /X2 controls the 
scaling behavior at Anderson transitions with short range 
interaction (ljj . In terminology of Wegner who pioneered 
the Anderson- localization multifractal analysis 17( , these 
exponents are denoted as X2 S < and x^ a > 0, respec- 
tively. To the linear order in couplings 7$, the RG equa- 
tions of the interacting theory should be controlled by 
scaling dimensions at the non-interacting fixed point. We 
have verified that this is indeed the case. Specifically, the 
exponent A = 2t, which is the positive eigenvalue of the 
matrix in ([1]) and which shows up in Eq. (O, is nothing 
but the anomalous fractal dimension (with an opposite 
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FIG. 2: (Color online) Temperature dependence of resistiv- 
ity t in the 2D symplectic symmetry class from numerical 
solution of Eqs. ©-© for 7^0 = -0.005, 7 c0 = -0.04, and 
t = 0.02,0.07,0.12,0.17 (from bottom to top). Left inset: 
projection of RG flow on a plane t = const. Separa- 

trix (thick black curve) approaches the line 7„ = 2y c . Right 
inset: dependence of T c on the bare resisitivity to- 



sign), A = — A2, while the second eigenvalue is the irrel- 
evant exponent, A' = — t = Thus, the enhancement 
of T c is intimately related to multifractality. 

Let us consider now a 2D system with strong spin-orbit 
interaction. In this case the spin-rotation symmetry is 
broken and the system belongs to the symplectic symme- 
try class. The change of the symmetry leads to two im- 
portant modifications of RG equations: (i) WAL replaces 
WL, and (ii) triplet interaction channel gets suppressed 
and can be discarded. The equations for the remaining 
interaction constants 7 S , "f c and resistivity t read 



dy V 7. 
dt/dy 
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The solution of Eq. (0 describes the WAL flow 
t- 1 (y)=t- 1 +y/2. 



(8) 
(9) 

(10) 



The eigenvalues of the linear part of the system (|8]) are 
again related to the multifractality of wave functions (this 
time for the symplectic symmetry class): A = t/2 = — A2; 
A' = — t = — fi2- The dominant eigenvalue A corresponds 
to the direction — 7 S = j c = 7, see Fig. [2] Projecting the 
system onto this eigenvector, we get 



d 1 /dy = (t/2) 7 - (4/3) 7 2 



(11) 



Assuming 70 = (—7^.0 + %lc,o)/3 < and solving 
Eq. (|TT|) . we find the scale at which the coupling (7! be- 
comes unity. This gives the transition temperature 



T, 



'Cxp{-C/(i |7o|) 1/2 } 



(12) 



where a numerical prcfactor C ~ 1 depends on the ratio 
7c.o/7s,o- Equation (|12[) is valid for I70I <C to; in the 



opposite case the clean BCS result (J?) is restored. The 
enhancement of T c in Eq. (fT2"]) is again due to multifrac- 
tality represented by the eigenvalue A = t/2 = — A2 in 
Eq. (jllj) . The enhancement is less efficient as compared 
to the orthogonal symmetry class, Eq. ©, because of an- 
tilocalizing behavior that leads to decrease of t and thus 
to weakening of multifractality (see Fig. [2]). 

We consider now a system at the Anderson transition 
point. This may be either a 2D or 3D symplectic class- 
system, or a 3D system of orthogonal symmetry. In all 
the cases, after the initial (fast) part of the RG evolution, 
where the Cooper term 7^ is unimportant, 7 S and (in 
the orthogonal case) 74 "adjust" to 7 C according to 7 S = 
— 7t = —7c (orthogonal) or j s = —~/ c (symplectic). So, 
the main part of the RG evolution can be described by a 
single equation, in analogy with Eqs. (|5|) and (|11[) : 



dj/dy = — A27 — a"f 2 



1 



(13) 



where A2 < is the fractal exponent for the given tran- 
sition point. In particular, A2 = —1.7 ± 0.05 for the 3D 
orthogonal-class and A 2 = -0.344 ± 0.004 for the 2D 
symplectic-class Anderson transitions [l6|. The SC will 
take place if 70 < 0. Analyzing Eq. (|13p . we find 



T* 



l7o 



|d/|A 2 | 



(14) 



for the transition temperature in d spatial dimensions. 
We see that at the Anderson transition point the en- 
hancement of the superconducting T c becomes even 
stronger than in 2D: T c is now a power-law (rather than 
exponential) function of the interaction constant. Equa- 
tion (fT4|) agrees with the result of Ref . [12] ■ 

The RG analysis allows us also to analyze the situa- 
tion when the system is slightly off the Anderson transi- 
tion. It is convenient to characterize the distance to the 
critical point i* by the correlation (localization) length 
£, ~ |*o — t* \~ v ■ The corresponding energy scale is the one- 
particle level spacing in the correlation volume 5^ oc l/£ d . 
The result fTJJ retains its validity in a vicinity of as 
long as <$£ <T C . On the insulating side (to > £*) the condi- 
tion (5j ~ T c determines the point of the superconductor- 
insulator quantum phase transition (see Fig. [3]) ■ On 
the metallic side (to < t*), there is a crossover regime 
extending from 6% ~ T* to 5% ~ 1 that provides a match- 
ing between the result fTJ]) at the Anderson-transition 
critical point and the clean BCS result 



14|. 



We close the paper with several comments: 

1) The true superconducting transition in 2D is of 
Berezinskii-Kostcrliz-Thouless (BKT) character, whereas 
we have calculated the mean-field transition temperature 
T c . It was found however [2l[ that the corresponding 
temperatures do not differ much, Tbkt — T c . Therefore, 
our result that shows an exponential enhancement of T c 
is expected to hold for Tbkt as well. 

2) The key assumption of the above theory was the 
neglect of long-range (1/r) Coulomb interaction. We 
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FIG. 3: (Color online) Schematic phase diagram in the inter- 
action-disorder plane near the critical point. Solid (blue) line 
denotes SIT. Dependence of T c on the distance from t* at fixed 
value of 70 is shown by solid (black) curve. In the quantum 
critical regime (region I) T c is given by Eq. (|14p . Away from 
criticality at t<t» (region III) the BCS expression @ holds. 
In the crossover regime (region II) T c ~£ -3 exp(— C3^ 2 /l7c.o|) 
in 3D and T c ~ ex P(~ c 2 a/? A2 / |Tc,o | ) for the symplectic 
symmetry class in 2D [Til ]. 



can think of the following situations when this should 
be justified: (i) a 3D material with a large background 
dielectric constant e; (ii) a 2D system on a susbstrate (or 
between two dielectrics) with large e; (hi) a 2D system 
with a nearby screening metallic layer; (iv) a system of 
interacting neutral fcrmions (i.e. cold atoms). 

3) A natural question is whether this effect has al- 
ready been observed in experiment. While one does see 
some enhancement of SC by disorder in several materials, 
one needs to argue that the Coulomb interaction is sup- 
pressed in order to attribute the increase of T c to the ef- 
fect of localization. In particular, a non-monotonous de- 
pendence of T c on the normal-state resistivity was found 
in Refs. [22, 23 1 in structures for which large e is expected. 

To summarize, we have developed a cr-model RG the- 
ory describing interplay of SC and Anderson localiza- 
tion in a disordered system with short-range interactions. 
This theory predicts a strong enhancement of SC by An- 
derson localization in 2D systems (at intermediate disor- 
der) and near localization transitions, implying a strongly 
non-monotonous dependence of T c on normal-state resis- 
tivity. Remarkably, the localization physics responsible 
for increase of resistivity and thus driving the system to- 
wards an insulating state favors at the same time the SC. 
It remains to be seen whether this mechanism may be 
employed in practice to obtain structures with strongly 
enhanced T c . The key condition is a suppression of the 
long-range component of the Coulomb interaction. 
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SUPPLEMENTARY INFORMATION 



I. ONE-LOOP RENORMALIZATION IN THE FINKEL'STEIN NLSM 



A. Nonlinear a-model: Definitions 



The action of the Non-Linear Sigma Model (NLSM) is given as a sum of the non- interacting part, S a , and contri- 
utions arising from the intcracl 
(Cooper), S$., channels [S1,S2]: 



butions arising from the interaction in the particle- hole singlet, S^, particle- hole triplet, S^ t , and particle-particle 



S — S& -\- >Sj 



(p) 



c( CT ) I QV 
int 1 °int "r~ a int , 



-Ar) 



where 



Sa = J drTr ( V( 3) 2 + 47rTz J drTr?](Q-A) 



,(p) 



°int 



32 

a, 7i r— 0,3 



drTr 



Tr 



a,n r— 0,3 j — 1 



T r t E E E drTl WrjQ Tr i a _ n t r jQ 



-Ac) 



a,n r=0,3 

Here g is the total Drude conductivity (in units e 2 /h and including spin) and we use the following matrices 
Ki = S gnn6 nm 5 a % , = n<W a/3 *oo, = 5 n - m , k 5^5^t m 

with a, /3 standing for replica indices and n, m corresponding to Matsubara fermionic energies e n = irT(2n + 1). 
The matrices 



t 



(15) 

(16) 
(17) 

(18) 
(19) 

(20) 

(21) 



operate in the particle-hole (index r) and spin (index j) spaces, with the corresponding Pauli matrices denoted by 

(22) 

(23) 



TO 



1 
1 



n 



so 



1 
1/ ' 



si 



r 

1 oj ' 

r 

1 o, 



T-2 



-i 

1 



T"3 



1 

-1 



S2 



— i 

1 



S3 



1 

-1 



The matrix field Q(r) (as well as the trace Tr) acts in the replica, Matsubara, spin, and particle-hole spaces. It obeys 
the following constraints: 



Q 2 = l, TrQ = 0, Q t = C T Q T C, 

where 

C = it 12 , C T = -C. 
It can be useful to represent Q as Q = T~ 1 AT, with matrices T obeying 

CT* = -TC, (T-^C = -CT' 1 . 



(24) 
(25) 
(26) 



In order to avoid notational confusion, it is instructive to compare our notation with that of the reviews by 
Finkel'stein [SI] and by Belitz and Kirkpatrick [S2]. First of all, these authors use different definitions of Pauli 
matrices: 



Rcf. [SI] 
Ref. [S2] 



F 

T a = T o, 

~BK 



TO, 



BK 



ITj, 



= So, (Tj=Sj, 



„BK 



SO, 



„BK 



-ISj 



j = 1,2,3 
J = 1,2, 3 



(27) 
(28) 
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The interaction terms ([T7|). (JTSJ) and ([19]) coincide with terms in Eqs. (3.9a), (3.9b), and (3.9b) in Rcf. [SI] with the 
coupling constants 

r. 5 = - T ^, r t = T r 2 , r c = T r c , (29) 

and with terms in Eqs. (3.92d), (3.92e), and (3.92f) in Ref. [S2] with 

r a = KW, r t = K<> 2 \ t c = k^/2, (30) 

where v is the thermodynamic density of states including spin. Finally, the parameters g and z in Eq. (|16[) are related 
by g = AttvD and z = (nv/A)Z to the corresponding parameters introduced in Ref. [SI] and by g = 1Q/G and z = H/2 
to those in Ref. [S2]. Note that Ref. [SI] focuses on the case of unscreened (long-ranged) Coulomb interaction: hence 
the interaction amplitude T s in the singlet particle-hole channel is expressed through the frequency renormalization 
factor Z there. In our case of a short-range interaction, these quantities are independent variables. 

B. Perturbation expansion 

We shall use the square-root parametrization 

q^w + aVi-w^ w ={ltj- ( 31 ) 

We adopt the following notations: W„ iri2 — w„ in2 and W„ 2ril = w„ 2ni with n% > and n 2 ^ 0. The blocks (in 
Matsubara space) obey 

w = -Cw T C, w = -Cw*C. (32) 

The second equality implies that in the expansion u>"f„ 2 = Y^rj( w nin 2 )rjtrj some of the elements (w"^ n2 ) are real 
and some are purely imaginary. 

Expanding the action S a to the second order in W, we find 



^ = ~\ E E E / jZl + [<J*[<CJ,, 03) 

rj ni»2 otf3 • ' 



/ d P 

Hence, the propagator becomes 



(K^lf.iifeU) = lD p (n 12 )S a ^S^S nin3 S n2n J rir2 S jlh , (34) 

D-\n 12 ) =p 2 + n 12 . 35 
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Note that (ww) and (ww) are also non-zero. From Eq. (|35[) we find 

(TrawTrbw) = y E E/ (^d^^ 12 ) E^nfni^n^lrj « 4F Tr -y-"-^— 6 ' 



(36) 



In what follows we will use the following identity: 

(Tr AW Tr BW) = 2Y Tr [AB - AAAB - ACB T C + AAACB T C\ . (38) 

C. Comments on interaction in the Cooper channel 

The Cooper channel interaction term S-£j. = (—ttTT c /2)0^ can be rewritten as 
1 3 

° (C) = o E E £>MSQ] ^LnQ] , {Ktl = S k+m>n 6^ 6^ t 00 (39) 



2 

a,n r— 1,2 j— 



However, Tr[t rj L%Q] = for j = 1, 2, 3 since 

Tr[t rj LZQ] = Tr[t r jL"Q] T = -Tr[Ct^CL a n Q] = -Tr[Cl* rj CL°Q] = J' 2 g ^[^Q] ' (40) 

Therefore, the operator 0^ c > describing the interaction in the Cooper channel is fully determined by the Cooper-singlet 
channel: 

° (C) = ^EE T r[^Q] Tr[t r0 LZQ] ■ (41) 

a,n r— 1,2 

D. Background field renormalization 

1. Singlet interaction in particle-hole channel: S^l = (— 7r2T s /4)0' p ' 
Performing transformation Q — > T^ 1 QTq, we find with the help of Eq. ([35]) and (|2"rJl) : 

° (P) = E E Tr [OroQ] Tr[/f„t r0 g] -^(E E ^[-^W^To] Tv[r n t rQ T^QT Q }\ 

a ? n r— 0,3 \ a,n r— 0,3 / 

= \E E Tr[T Z n Q ^oT ( 7V]Tr[To/f n i I .or o - 1 ^]\ =2F^ E Tr[-I n % Qo/- r AoQo + I%t r0 Q I%Ct* r0 CQ ] 

\ a,n r— 0,3 / a,n r— 0,3 

= - 2F E E {Tr[/:t r oQor„t r0 Qo] + Tr[lZt r oQ i:t r0 Qo\}. (42) 

a,n r— 0,3 

Next, we represent the resulting expressions in terms of the operators entering in the interaction- induced part of the 
action, ([17]), (HHJ), and dHJ): 

E E Tr[l%,t r0 QoI« n .t r oQ ] = \ E E E ^[WnQo] Tr[l« n t rj Q ] = + (43) 

a,n' r— 0,3 ct,n r— 0,3 j—0 

Finally, we see that the renormalized particle-hole singlet term is expressed through all three interaction operators : 

O(p) 0(p) _ YO {p) - YO {a) - 2YO {c) . (44) 

2. Triplet interaction in particle-hole channel: S^2 = (— 7rTTt/4)0^ 
Performing transformation Q — > T^ 1 QTq, we find with the help of Eq. (|38|) and (|26j) : 

= E E ETr[/«^g] Tr^t^Q] ^ E E^W^o] Tr [r n t rj T^QT ] \ 

Q,n r-0,3 j-1 \ a.n r— 0,3 j — 1 / 

= (E E Y.H T ^nU ] T^W]TT[T^ n tr 3 T^ 1 W]\ 

\ a,n r— 0,3 j — 1 / 
3 

= 21^E E E Tr[-/"t rJ Q -ff n^rjQo + In^rjQaln Ct*jCQ ] 
a,n r— 0,3 J — 1 
3 

= - 2Y E E E{^[ 7 ™*^° J -«^o] - (-If Tr[0,,g O rj Qo] }■ (45) 

a,n r— 0,3 j—1 
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Next, 



E E ^^[I^UjQoI^trjQo] = ^E E {3Tr[O r0 Q ] Tr[I« n t r0 Q Q ]-J2Tr[l%t rj Q Q ] Tr[I« n t rj Q ] 

a,n' r— 0,3 j — 1 a,n r-0,3 j — 1 

= \0 {p) -\0 { °\ (46) 



and 



E E (-l) r ETr[Z«i rj Qo4 Q ,^Oo] = ^E E 3Tr[L^ r0 Q ] Tr[L«i r0 Qo] - ^ Tr[L«t r ,Q ] Tr[i^Q ] I 

a,n' r— 0,3 j — 1 a.n r—1,2 \ j — 1 ) 

= I E E Tr [ L »^oQo] Tr[L"t r0 Qo] = 30 (c) . (47) 

a,n r— 1,2 

Finally, we find 

i") -j. (ct) - 3FO (p) + rO (CT) + 6YO (c) . (48) 

5. Interaction in Cooper channel: S$. = (-7r2T c /2)0 (c) 
Performing transformation Q — > Tq^QTo, we find with the help of Eq. (|38p and (f2"cJl) : 

° (C) = E E (-l) r Tr[/«t r0 Q/« ri i r0 g] = ^E E ^[^*roQ] Tr[LZt r0 Q] 

a,n r— 0,3 a,7i r— 1,2 

^ (E E ^[iSW'QTo] Tr^^oTQ-^To]) 

\ a,n r=l,2 / 

J(E E TrfToi^oTo-^jTrfToL^oTo-V] ) =y^ ^ Tr[-L,% Qo£^oQo + L^ r oQo^Ct; o CQo] 



2 



l a.n r— 1,2 / a,n r— 1,2 



= -2YYl E ^[OroQo^troQo]- (49) 
a,n r— 1,2 

Next, 

1 3 

53 ^ Tr[L«^Q L^,oQo] = ^E E (-l)'X Tr [^Qo] Tr[O rj Q ] 

a, r;/ r— 1,2 a,n r— 0,3 j—0 

i 3 

= "2 E E (-l) r E Tl 10^o] Tr[/« n C^.Cg ] 

a,n r— 0,3 j—0 

= -jEE j ^[OroQo] Tr[l«J r0 Q ] - ^ Tr[O r ,Q ] Tr[r n t rj Q ] ] = -io<"> + . (50) 

q,ti r=0,3 I j=l ) 

Finally, we obtain 

0(c) ^ 0(c) _ y (p) + YO {a) . (51) 

It is worthwhile to mention that only particle-hole singlet and triplet interaction operators are involved in the renor- 
malized particle-particle term. It is due to the different properties of and L" matrices under the transposition 
operation: (7«) T = (L«) T = L° 
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E. One-loop renormalization group equations 



Using Eqs. (|34"]). (|48)l and (|5T|), we find 



Next, since 



~(r t - T S Y + T t Y + 2T C Y)0^ - Z^(r c - T s Y + 3T t Y)0^ (52) 



g J (2TT) d irg 



we obtain the one-loop equations for the renormalization of the three interaction coupling constants 

= -— (r s + 3r t + 2r c ), (54) 



d In L irg 

#V = -— (r.-r t -2r c ), (55) 

a In L Trg 

#y = --(r s -3r t ). (56) 

a In L -Kg 

The right-hand sides of these equations are linear in the interaction couplings, since we have neglected the ladder 
resummation in the interaction propagators, as appropriate for weak couplings. The flow of the field renormalization 
constant z follows from the renormalization of T s in view of the particle number conservation dT s /dhiL = —dz/dhiL: 

dz dT, 1 , , . . 

Introducing "f s ,t.c = ^s,t,c/z and taking into account the definition of t = 2/ng, we arrive at the linearized RG 
equations of the main text [see equations (|70)) -([T4" )) below]: 

^( 7s + 3 7t + 2 7c ), (58) 



din I, 2 

= -kls-lt-2 7c ), (59) 



dlnL 2 

k-Y.-*rt)- (60) 



dlnL 2 
F. Relation to the BCS Hamiltonian 



Let us consider the interaction part of the Hamiltonian: 



H int = i / drdr'U(r - r')Vv WVvMVv (OVv (O (61) 



2 

In the BCS case (for example, for short-range attraction mediated by phonons), 

U(r - r') = --S(r - r'), 
v 

where the thermodynamic density of states v accounts for spin. According to Ref. [S2] the interaction parameters can 
be written as 

7, = -itW 7s = -lffe' lc = ' Fc (62) 

where 

F t = ~(U(2k F sm(6/2))) FS , F s = vU(q) + F t , F c = ~{U{2k F sm{d/2))) FS - ^(U(2k F co S (6/2))) FS . 

(63) 
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Here {■ ■ - }fs denotes averaging over the Fermi surface. In the case of BSC Hamiltonian we find 

F t = A/2, F s = -A/2, F c = A/2 (64) 

and 

7s«A/2, 7c = -A/2. (65) 

Thus, for the BCS case (i.e. when neither screened nor unscreened Coulomb repulsion is taken into account), we 
get the following interaction parameters at the ultraviolet scale (which is given by Dcbye frequency ujp in the case of 
phonon-induced superconductivity) : 

-7s = It = 7c = -A/2. 

As we will see in Sec. II below, this is precisely the relevant direction for the RG flow. When disorder is weak, 
uipT 3> 1, the initial values for these couplings in Eqs. ([58 j) -([60" j) are taken at the scale corresponding to the elastic 
scattering rate 1/t. Then the Cooper interaction constant is renormalized at ballistic scales (between 1/r and wj>) 
such that 

-^o=7t,o = -A/2, Jcfi = - —^-^ = -A uw - (66) 

where T C BCS = lo d exp(-2/A). 

The RG equations imply that there is the following pcrturbative correction to the BCS coupling due to interaction 
in the particle-hole channel 

7c = 7c,o - -Y^JTi ( 67 ) 
which is logarithmic in temperature. This leads to increase of the transition temperature: 

ST C _ t Q X 1 _ t X ( 1 V ( . 

This perturbative calculation is justified in the regime 5T C /T® CS <C 1 which corresponds to the condition to <C 
l7c,o| 3 /A. In order to neglect the renormalization of t due to weak localization disorder should be sufficiently weak, 
to <C |7co (see Sec. II). Therefore, the perturbative correction to T c given by Eq. (155)) is valid for to/|7co| ^ 
min{l,| 7 ; | 2 /A}. 

In sufficiently dirty systems, lue>t <C 1, the bare values of interactions constans in Eqs. (|58l) - (|60l are taken at the 
scale corresponding to wd. In this case, 

ST C 2t Q u D 4i /cn s 



TBCS X TBCS x 2 

which holds for to ^ 7c o- I n * ne next section we will analyze the case of stronger disorder for which the perturbative 
treatment is insufficient and one should solve the full set of RG equations. 



II. ANALYSIS OF RENORMALIZATION GROUP EQUATIONS 



We consider first 2D systems of the orthogonal symmetry class at weak disorder (g ^> 1) where the localization and 
the multifractality effects are weak. Then we turn to 2D systems of the symplectic symmetry class in the metallic 
regime (g 3> 1) characterized by weak antilocalization. Finally, we address the vicinity of an Anderson transition 
(symplectic class in 2D and orthogonal class in 3D). In all these cases the multifractal character of wave functions 
in dirty systems can strongly enhance the superconducting transition temperature as compared to that of the clean 
system (usual BCS). 
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A. 2D, orthogonal symmetry class 

The full set of RG equations for the orthogonal symmetry class reads: 

dt ,2 (Is , „7t , \j 



-^(7 s +3 7t + 27c) I (71) 



dlnZ 2 V 

~(7.-7*-27c), (72) 



din I, 2 

(7 s -3 7t )-2 7c 2 , 



d 7c f,_ 



dlni 2 
dlnz t 



d\nL 2 



(7. + 37* + 2 7c ) • (74) 



Here, we remind, t <C 1 is the dimcnsionlcss resistance (inverse dimensionlcss conductance, t = 2/irg), ji = Yi/z <C 1 
are interaction constants in the singlet (s), triplet (t) and Cooper (c) channels. In Eq. (|70| the first term is due 
to weak localization. The second term is due to interaction [Altshulcr-Aronov (AA) correction], it can be neglected 
within our accuracy. In Eqs. (|71 |) -(|73 j) only terms linear in interaction constants have been kept. The only exception 
is the Cooper term — 2 7 ^ in Eqs. (|73|) which is responsible for the BCS superconducting transition. The neglected 
terms are much smaller because they contain an additional small factor t. Below we shall use y = In L where L is the 
running RG length scale. We measure lengths in units of the microscopic scale where the RG ([70]) - ([74)) starts. 

Neglecting the AA contributions, we rewrite the set of equations (fT0)) - (|73|) in the form 



(75) 

z V 1 -3 J \ 7c / \ 2 7c 2 / 

Ty =t - (76) 

Eigenvalues of the matrix in Eq. (|75|) arc A = —4, A' = 2,2 (not including the prefactor — i/2); the corresponding 
eigenvectors are as follows: 



A = -4 : 1 ; A' = 2 : 1 and -1 . (77) 








If the 7^ term is neglected, the solution of the linear system (|75[l approaches the eigenvector with A = —4, i.e., 
7s — —"ft = — 7c- Let us expand the vector formed by 7 j in eigenvectors (|77|) : 



(78) 




-1 


1 






1 


1 




• 


1 


-1 







For convenience, we also present the inverted transformation: 

-1/6 1/2 1/3 

1/2 1/2 

1/3 1/3 







• 


H 








(79) 
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Transforming the set of equations (|75[) to the new variables a, b, c, we get 

da „ 2 



2ta- =-(a-b + 2c) 2 , (80) 
ay 3 



= -tb, (81) 



d6 

- = -tc--(a-b + 2c) 2 , (82) 



dt 
dy 



r . (83) 



Equations ([80| -(|83 |l are supplemented by the following initial conditions: a(0) = ao, 6(0) = M c(0) = Co and i(0) = to- 
Equation ([83]) yield the standard evolution of the resistance due to weak localization: 



Equation ([8lj) can now be solved, yielding 



t(y) = ■ (84) 
i - toy 



b(y) = Ml - toy) = 6 y • (85) 



Since b decreases upon RG, it is not important and we neglect it in the future analysis. 

Equations for the remaining two variables, a and c are coupled. If the quadratic term is neglected, then a increases 
and c decreases. This suggests that c can be neglected. This is confirmed by a more careful analysis which shows 
that, although on the very last interval of RG "time" y the variable c starts to increase and becomes of the same 
order as a (i.e. of order unity), this weakly affect the RG scale at which this happens (i.e. the temperature of the 
superconducting transition). Thus, we neglect c in what follows. 

We can now easily solve the remaining equation for a. We assume the starting value ao to be negative (which means 
that there is attraction in the Cooper channel that is supposed to lead to the superconductivity), a = —\a\. This is 
in particular the case when 7 Ci o is the dominant coupling and 7 Ci o < 0. The equation reads 



M = 2 t\a\+ 2 -a\ (86) 
dy 3 



Solving this equation, we obtain 



Let us analyze the obtained result. Let us first assume that |ao| >C to- Then the second term in brackets in the 
r.h.s. of (|87|l is small compared to the first one and can be neglected, 

2 



yy > t\\ao\t 3 J 



With increasing RG scale y the resistance t increases together with the interaction a. If t reaches first unity, we 
get an insulator; if a ~ 1 happens first, we get a superconductor. It is easy to see that the second possibility 
(superconductivity) is realized if |oq| 3> t\. Then at the point of the transition to superconductivity we have a 
resistance 

% « 1 . (89) 



2|ao 



The transition will then happen at 
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i.e. at temperature 

(91) 

Here the factor 2 in the exponent originates from a translation of the length scale into energy (temperature) . Under 
the above assumption |oq| <C to the second term in square brackets in the exponential of (|9ip is just a small correction 
to the first one. By solving (|82|) with b = and a given by Eq. ([88]) . one finds that although |c| decreases initially, 
eventually with increasing RG scale towards y* it becomes of the order of unity: |c(j/*)| ~ 1. Therefore, to determine 
precise value of t* one has to solve coupled equations for a and c (with (b = 0)). 

The transition temperature (|9ip is much higher than the BCS temperature Tbcs ~ e" 1 /'^ ', so that the super- 
conductivity is strongly enhanced by disorder. The origin of the enhanced superconductivity is in the increase of \a\ 
governed by the eigenvalue A = — 4 of Eq. ([77)) which yields the eigenvalue — (i/2) x (—4) = 2t of the linear part of 
the system (|75[) . This is nothing but the anomalous multifractal exponent — A 2 for this symmetry class. (We have 
in mind the "weak multifractality" in 2D.) Therefore, the (multi)fractality is the source of the enhancement of the 
supcroenductivity. 

It is worth mentioning that there is no such enhancement in the absence of interaction in particle-hole singlet 
and triplet channels. Indeed, the interplay of disorder and interaction in the renormalization of j c does not produce 
the term t^ c on the r.h.s. Eq. (|73|) . Of course, if initially absent, the triplet and singlet amplitudes are generated 
due to the interplay of the Cooper amplitude and disorder, see Eq. (fTTj) and (|72|) . Furthermore, as we have seen in 
Sec. IF above, the BCS interaction in fact contributes to all interaction channels, so that even in the absence of 
electron-electron repulsion the bare values of 7 S and 7t are non-zero. 

If |ao| <C ig; the resistance reaches unity before the interaction becomes strong, and the system is an insulator. 
Finally, if |c*o | 3> to, the disorder is not particularly important, and the transition temperature is given by usual clean 
BCS, T t ~ e-VWc.ol. (i n tne latter case neglecting b and c is not paramctrically justified and leads to an incorrect 
numerical factor in the exponent.) 



t;~cx P {-A i 



B. 2D, symplectic symmetry class, metallic regime 



Now we consider the symplectic class, i.e. assume that the spin symmetry is broken for example, by sufficiently 
strong spin-orbit interaction. This leads to the following two modifications: (i) weak antilocalization rather than weak 
localization, and (ii) triplet interaction channel is suppressed and can be discarded. With these modifications, the 
system (|7Q]h(|73|) becomes (cf. Ref. [S4]) 



dt 
dy 



-r = -o* - Hf+7c t 



d'fs 
dy 

djc 
dy 



-2(7. + 27c), 



(92) 
(93) 
(94) 



Similar equations for the symplectic case of long-range Coulomb interaction have been derived in Ref. [S5]. The 
Altshuler-Aronov terms in (ji?2"j) can again be neglected. Equation (|9"2"|) then yields the standard antilocalization 
behavior, 



t(y) 



to 



1 + yt /2 



The system of equations (|92|) . (|93|) in the matrix form is: 



dy V 7c 

Eigenvalues of the matrix in Eq. (|96|) are A 
eigenvectors are as follows: 

A = -1 : 



1 2 
1 





27c 2 



(95) 



(96) 



-1, A' = 2 (not including the prefactor —t/2); the corresponding 



A' = 2 



(97) 
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(98) 



(99) 



If the 7^ term is neglected, the solution of the linear system approaches the eigenvector with A = — 1, i.e., 7 S 
As in the orthogonal case, we can expand the vector formed by 7 j in eigenvectors 

The inverted transformation is given as: 

-1/3 2/3 \ / 7s 
1/3 1/3 J { 7c 

Transforming the set of equations (|98[) to the new variables a, c, we find 

^=-tc-\{a + cf. (101) 
ay 3 

Equations for two variables, a and c are coupled. If the quadratic term is neglected, then a increases and c decreases. 
At the later stage of RG the quadratic terms leads to enhancement of c. This suggests that c can be neglected for 
qualitative analysis of RG equations (|100|) - (|101[) . We thus neglect c and keep only a (fully analogously to what we 
have done in the orthogonal case). The resuting equation for a reads 



d\a\ t , 4 o , 

V = 2 lal+ r- (102) 



+ 4-4) 



We solve this equation with the result 

t \a t ' 3t§ 3t 2 

The new (different from usual BCS) behavior emerges under the condition ao "C to- Then the condition a ~ 1 yields 

U ~ 2 (a t /3) 1/2 < t Q . (104) 

By solving Eq. (|101[) with a given by Eq. (|103[) . we find that although |c decreases initially, eventually with increasing 
RG scale it reaches a: c ~ a at t = f». Therefore, to determine precise value of f» one has to solve coupled equations 
for a and c. 

Equation (|104[) yields the transition temperature 

T* ~ e~ 2 ^ ~ e- 4 /** ~ exp {-C/ Mo) 1/2 } , (105) 

much higher than Tbcs ~ e~ 1 ''" /c ^. The BCS behavior is restored (up to corrections) at ao ^> to. The constant C is 
of the order of unity and depend on ration co/ag- 

As in the orthogonal case, the source of the enhancement of the superconducting temperature is in the first term on 
the r.h.s. of eq. (|102|) . The eigenvalue t/2 is the anomalous multifractal exponent — A2 for the symplectic symmetry 
class. Therefore, also in this case the (multi)fractality is the source of the enhancement of the superoenductivity. This 
enhancement is less efficient than in the orthogonal case for two reasons, because of antilocalizing behavior that leads 
to decrease of t and therefore weakening of multifractality. 



C. System at or near Anderson transition 

1. Exactly at criticality 

We now consider a system at the Anderson transition point. This may be 2D or 3D symplectic class system, or 3D 
orthogonal class. In all the cases after the first (fast) part of the RG evolution, where the Cooper term is assumed 
to be not important yet, the coupling constant 7^ and (in the orthogonal case) 7 t "adjust" to 7 C according to the 
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7s = — 7< = — 7c (orthogonal) or 7 S = — 7 C (symplcctic). So, the main part of evolution can be described by a single 
equation, in analogy with Eqs. (f86|) and f|102[) : 

P- = -A 27 - a 7 2 , a - 1 (106) 
ay 

The superconductivity will take place if 70 < 0. We also note that A2 < 0. In particular, A2 = — 1.7 ± 0.05 for the 
3D orthogonal-class [S6] and A2 = — 0.344 ± 0.004 for the 2D symplcctic-class Anderson transitions [S7]. If I70I <C 1, 
the second term in the r.h.s. of (|106[) is not important for the evaluation of the leading behavior of the transition 
temperature. Keeping only the first term, we find 

T* ~ «p{-dy,} ~ | 7 o| d/|A2 ' • (107) 

The factor d in the exponential, where d stands for the spatial dimensionality, translates length in temperature. We 
see that at the Anderson transition point the enhancement of the superconducting transition temperature becomes 
even stronger than in 2D: the transition temperature is now a power-law (rather than exponential) function of the 
coupling constant 70. Equation (|107[) agrees with the result of Rcf . [S3]. 



2. Slightly off criticality: Insulating side 

In the above consideration we assumed that the system is exactly at the Anderson transition point. Let us analyze 
what happens if the system is slightly off criticality. Then, we need to add the following equation to Eq. (| 106(1 : 

dt 1 , , 
— = -{t-t c )+ m . 108 
dy v 

Here, we take in account that the presence of interaction drives the system away from the non-interacting critical 
point. We note that the correlation length exponent v > 0. In particular, v = 1.57 ± 0.02 for the 3D orthogonal-class 
[S8] and v — 2.746 ± 0.009 for the 2D symplectic-class Anderson transitions [S9]. 
Neglecting the term of the second order in 7, the system of Eqs. (|106|) and (|108[) can be solved for |A 2 |^ 7^ 1 as 

t = t+ Y 7, t = t c +(t -t c )e^, 7 = 7oe l A ^, t =t - -^ (109) 
|A 2 |j/ — 1 V J \J\ 2 \v — 1 

In the special case \&2\v — 1, we find 

t = t-r)j jy, t =t . (110) 

Therefore, the presence of the term 777 in Eq. (|108p indicates that the proper scaling variable is t rather than t. In 
what follows we shall omit 'tilde' sign. 
Up to the scale of the localization length 

Z~\to-t c \- v (111) 

the RG will proceed as at the critical point. So, there are two possibilities. If L* < £, where L* ~ 7o|~ 1// ' A2 
is the length scale where the superconducting transition at criticality takes place (Sec. II CI), then the transition 
temperature is not affected by detuning. In the opposite case, the localization takes place first, and there is no 
superconductivity. Thus, the condition (6% oc £~ d ) 

5 e ~ T* , (112) 

where T* is given by Eq. (|107[) , is the condition of the superconductor- insulator transition (at zero temperature). It 
is worth noting that Ref. [S3] argues that supeconducting state with T c <C T* persists further in the localized regime 
(<5j ^> T*) due to Mott-type rare configurations. Our RG approach (at least in its present form) is not sufficient to 
explore this posibility. 



3. Slightly off criticality: Metallic side 

As in Sec. II C 2, we have the length scale £ given by Eq. (jllip but it now has a meaning of the correlation length. 
Below this scale the system is at criticality, above this scale it becomes a metal. As on the localized side, if <5j «!„ 
the detuning from criticality is immaterial, and the transition temperature is given by T*, Eq. (|107[) . In the opposite 
case, the result depend on whether we are in 2D symplcctic case or in 3D (4D, ...) orthogonal class. 
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4- 3D, 4D, . . . , orthogonal class 

After the first step of RG (up to the scale £) the Cooper-channel interaction constant takes the value 

7o=7o£ |A2 '. (113) 

After this the RG proceeds according to the usual BCS. The total RG scale y where the coupling becomes of order 
unity (and thus the transition happens) is 

^=^+]^iSP ( 114 ) 



where c ~ l.Thus, the transition temperature is 



cd 



|A 2 |AT 



T - a ~r^\-¥^] =St ' xp {- eJ {f.) }" (115) 

It is easy to see that this intermediate regime correctly matches results for two regimes between which it is located: 
at ^ - T» we have tP ~ T*, and at ~ 1 we have ~ T BC s- 



5. 2D, symplectic class 

After the first (Anderson-critical) step of RG the interaction constant is given by the same formula (|113|) as in 3D. 
The difference is on the second step: now we have to apply the formula for a 2D symplectcic metal (|105l) with to ~ 1. 
This yields 

^~?^{-ra4rai}- 4, - p {- c (i)' 4 "'' 4 } ■ (116) 
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